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Abstract 

Geometrically, zeroes of a Gaussian analytic function are intersec- 
tion points of an analytic curve in a Hilbert space with a randomly 
chosen hyperplane. Mathematical physics provides another interpre- 
tation as a gas of interacting particles. In the last decade, these in- 
terpretations influenced progress in understanding statistical patterns 
in the zeroes of Gaussian analytic functions, and led to the discovery 
of canonical models with invariant zero distribution. We shall discuss 
some of recent results in this area and mention several open questions. 



Introduction 

A Gaussian analytic function is a linear combination 

of analytic functions fk : G — > C (G C C is a domain), 

\fk(z)\ 2 < oo locally uniformly in G, 

fc>0 

with independent standard complex Gaussian random coefficients Q^. The 
random zero set Zf = / _1 (0) is the theme of this talk. 

* Supported by the Israel Science Foundation of the Israel Academy of Sciences and 
Humanities. 



Pioneering contributions in this area were made by Paley and Wiener |28l 
Chapter X], Kac |151 Chapter I], and Rice [30. • Paley and Wiener con- 
structed a large class of Gaussian analytic functions in a strip with station- 
ary distribution with respect to the shifts, and computed the mean number 
of zeroes. Their work was influenced by ergodic theory and theory of almost- 
periodic functions. Kac was interested in the mean number of real zeroes 
of polynomials with real coefficients. Rice systematically treated both theo- 
retical and applied aspects of random noises in radio signals. The technique 
introduced by Kac and Rice has been significant to radio engineers and 
physicists. 

These studies were continued in various directions, notably by Little- 
wood and Offord [2Uj , Hammersley ^2] i Offord |261 127j , and Kahane |16l 
Chapter 13]. Hammersley looked at pure probabilistic aspects of point pro- 
cesses generated by zeroes of random polynomials. The other authors were 
motivated by the entire functions and Nevanlinna theory. Introducing ran- 
domness, they tried to single out 'typical patterns' in the zero distribution 
in various classes of entire and analytic functions. 

In the 90s, the subject was revived by several groups of researchers who 
came from different areas: Bogomolny, Bohigas and Leboeuf; Shub and 
Smale; Edelman and Kostlan; Hannay; Bleher, Shiffman and Zelditch; Non- 
nenmacher and Voros; by no means is this list complete. They established 
new links to physics (Coulomb gas of charged particles, random matrices, 
quantum chaos) and geometry (analytic curves in projective Hilbert spaces), 
and drastically changed the whole subject. 

We start this lecture with a quick review of basic results pertaining to 
zero sets of arbitrary Gaussian analytic functions (Section In Sectional 
we introduce three canonical random zero process (on the Riemann sphere, 
complex plane, and the unit disc), distinguished by stationarity with respect 
to the isometries. In Sections |H1 and 0J we consider in more detail one of 
them, the canonical random zero process in C. The exposition in this part 
is based on joint works with Tsirelson [32] • In Section we discuss an 
'exactly solvable case' of the hyperbolic zero process recently discovered by 
Peres and Virag |29j . 

1 The random zero process Zf 

Informally speaking, the random zero set Zf = / _1 (0) is the intersection of 
an analytic curve f : G — ► P{H) in a projective Hilbert space with a random 
hyperplane, the analytic functions fk are homogeneous coordinates of the 
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curve f. Projective unitary transformations of the curve f do not change the 
random zero set Z$. Hence the random set Z$ depends only on geometry of 
the curve f. Its study can be interpreted as part of the H. Cartan - Ahlfors 
- H. and J. Weyl theory of analytic curves independent of the dimension of 
the target space. 1 

The properties of the random process Zj can be expressed by its counting 
measure 

a: f(a)=0 

5 a is a point measure at a. The measure rif is a random, positive, locally 
finite measure on G. The classical formula 

n / = i-Alog|/| (1.1) 

(the Laplacian is understood in the sense of distributions) explicitly relates 
the random measure rif to the Gaussian analytic process /. Proofs of most 
of the results presented below start with this relation. 

1.1 The Edelman-Kostlan formula for the mean measure 

The first question about the random measure rif is to find its average which 
is a non-negative measure in G. 

Theorem 1.2 (Edelman-Kostlan |9j) 

En / = ^-Alog||f||, 

where 



y A:>0 

I.e., the mean measure Enj coincides with the Riesz measure of the subhar- 
monic function log ||f||(z). 

The RHS of the Edelman-Kostlan formula is a pull-back of the Fubini- 
Study area measure from the projective space P(H) to G by the curve f. Its 
density with respect to the Euclidean area measure in G equals 

1 (f#\2 — 1 ^i<fc \fjfk ~ fifk\ 2 (-> o\ 



An interesting attempt to build a 'dimensionless theory' of analytic curves was made 
by Favorov His approach is based on the pluripotential theory in Banach spaces. 
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The function is a 'Fubini-Study derivative' of the curve f. The Edelman- 
Kostlan formula can be viewed as a version of the classical Crofton formula 
from the integral geometry. Its proof is a simple computation based on 
equation Ql.lj) : 

En/ = i-A(E log |/|) = i-Alog H/ll, 

since, for any complex Gaussian random variable £, 

E log |C| = log HCII + const. 

What about the higher 'moments' of the random measure n/? They are 
expressed by the /c-point correlation measures 

dfj,(zi, ... z k ) = E (drif(zi) ... dn f {z k )) (1.4) 

on G x •• G - Hannay |13j derived explicit formulas for these measures which 

k times 

generalize (|1.3j) . They involve determinants and permanents of k x k matri- 
ces. In different contexts, the rigorous proof of these formulas is given in [3J 
and I2S1- 

1.2 Calabi's rigidity 

Surprisingly, the mean En/ determines the random zero set Zf. In geometry, 
the same phenomenon was discovered by Calabi already in the beginning of 
the 50s. 

Theorem 1.5 Let f and g be Gaussian analytic functions in a domain G, 
and let En/ = En 9 . Then the corresponding random zero sets Zf and Z g 
have the same distribution. 

This holds due to the underlying analyticity. The idea is not difficult: let 
K(z\,Z2) = E ^/ (z\) f (z2)j be the covariance of the process /. By the 
Edehnan-Kostlan formula, the mean measure En/ determines the function 
z i — ^ log K(z, z) up to a harmonic summand. In turn, the diagonal K(z, z) 
determines the whole covariance kernel K (due to analyticity of K in z\ and 
Z2), and hence the distribution of the Gaussian process /. The details and 
references are in |33j . 

Here is Calabi's original formulation: // two linearly non- degenerate an- 
alytic curves f : M — > P{H\), q: M — > P{H2) of a complex manifold M in- 
duce the same Riemannian metric on M by pulling back the corresponding 
Fubini-Study metrics, then the projective spaces coincide P{H\) = P(H2), 
and the curves are unitarily equivalent. 
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1.3 Offord-type estimate 



Theorem 1.6 Let f be a Gaussian analytic function on a plane domain G. 
Then for any test function <fi € Cq(G) with a compact support in G, and any 
A > 



a 



4>(dnf — E(dn/) ) 



> A I < 3 exp 



2vrA 



(1.7) 



Here, || . ||i is the L 1 norm with respect to the area measure. 

Here is an argument borrowed from Offord [2Sj- By (|1.1[) and Green's 
formula, we need to estimate the probability 



(log |/| -E log |/|) A0 dm 



> 2n\ 



m is the area measure. This reduces the proof to a simple fact about con- 
centration of log |C|, where £ is a complex Gaussian random variable. The 
details are in [3*3] . 

The result can be extended in various directions. It persists for zero sets 
of any random analytic process / in G with uniformly bounded exponential 
moment: for some c > 0, 



sup E ( e 

zeG 



c| lag l/toll 



< oo . 



Examples of such analytic processes are given in |24j . Instead of the L 1 - 
norm of the Laplacian A<p, one can fix any of the L 9 -norms, 2 < q < oo, 
of the gradient In Section |1J we discuss a more complicated 'global 

version' of Theorem 11.61 

There is a price for such a level of generality: sometimes, Offord's esti- 
mate does not give an optimal result. For example, it does not yield sharp 
bounds for the 'hole probability' (see (|3.4j) and ()5.3|) below). 



2 Chaotic analytic zero points 

In the beginning of the nineties, Bogomolny, Bohigas and Leboeuf; Kostlan; 
and Shub and Smale introduced a remarkably unique class of Gaussian an- 
alytic functions with unitary invariance of zero points. Following Hannay 
|13j . we use the term 'chaotic analytic zero points' (CAZP, for short). We 
consider here three CAZP models: the elliptic CAZP, the flat CAZP, and 
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the hyperbolic CAZP 2 . They are the random zero set of a Gaussian analytic 
function 




L(L - 1) . . . (L - k + 1) 
k\ 



■z k (elliptic, L = 1,2,...), (2.1) 
(fiat, L > 0), (2.2) 
h(z) = ^\l L{L + 1) " k < ; L + k ~ 1) * k (hyperbolic, L > 0). (2.3) 

fc=0 

The analytic function (|2.1|) is a polynomial of degree L (the domain of the 
elliptic CAZP is the Riemann sphere), the function (|2.2|) with probability 
one is an entire function, and the function (|2.3|) with probability one is 
analytic in the unit disc. 

We introduce unified notation: M for the domain of the CAZP, and 
r for the symmetry group of M. Then CAZP is a unique T -stationary 
random zero process. Here, T-stationarity means that for any 7 £ T the 
point processes Zf and Zfo-y have the same distribution. Uniqueness means 
that CAZP is the only T-stationary process on A4 among the random zero 
sets of Gaussian analytic functions. 

Having explicit formulas (|2.1j) . (|2.2j) . (|2.3j) . it is very easy to prove T- 
stationarity and uniqueness. It suffices only to check that 

En fL = L-m* , (2.4) 

where m* is a normalized T-invariant area measure on M. Then, by Calabi's 
rigidity, Zf L is T-stationary, and, again by Calabi's rigidity, fi is unique. 
Verification of (|2.4|) is a straightforward application of the Edelman-Kostlan 
formula. For example, in the flat case, ||/i|| 2 (z) = exp(L|z| 2 ), and 

n fL = ^A(L\z\ 2 /2)=L.-m. 

We see that in the flat case the normalized area m* equals —m. It is also 
worth mentioning that canonical isometric embeddings of M into projective 
Hilbert spaces corresponding to the Gaussian analytic functions ()2.1|) . (|2.2I) 
and (|2.3|) are well known in geometry and physics. 

2 The toric CAZP (which we do not discuss here) was introduced by Nonnenmacher 
and Voros |25j . 
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By (|2.4|) . the parameter L equals the mean number of random zeroes per 
unit area on Ai. In what follows, we shall consider the asymptotic behaviour 
of the random zero processes Zf L in the 'large intensity limit L — > oo. Some 
features for all three canonical models are similar, some are different. Due 
to compactness, the elliptic model sometimes is simpler to analyze. On the 
other hand, the hyperbolic model has additional intriguing features. 

To fix ideas, we shall concentrate on the flat model. In this case, intro- 
ducing L, we just make a homothety of the plane with coefficient r = \fL. 
This makes the flat case more transparent 3 . Thus, we do not need param- 
eter L anymore, and we consider the asymptotic zero distribution of the 
Gaussian entire function of order two 

OO 

k=0 

The function f(z) can be viewed as a Gaussian counterpart of the Weier- 
strass o"-function 

a(z) =z H (l - -) e*l«+m*lu? . 
LoeV^z 2 \{o} 10 

Indeed, the random function |/(z)|e~' 2 ' 2//2 has a stationary distribution, 
while the function |<r(2:)|e _ ' z ' ' 2 has periods y/rr and iy/n. 



3 Linear statistics 



Given a test-function h: C — > M with a compact support, consider the ran- 
dom variable 



Z r (h) = J h(l) dn(z) , EZ r (h) = r - J 



2 

hdm , 



n is a counting measure of the flat CAZP with intensity L = 1, m is the 
area measure. We are interested in the asymptotic behaviour of Z r (h) when 
r — ► oo. The size of fluctuations of Z r {h) depends on the smoothness of the 
test-function h. 



3 The scaling z — zo + flattens out the elliptic and hyperbolic geometry as L — » oo. 
In this limit, the entire function 12.21 is a locally uniform limit of the functions 12. It and 
and the flat CAZP appears as a scaling limit of the other two CAZP models. This 
is the motivation for an advanced theory developed by Bleher, Shiffman and Zelditch in 

m 
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3.1 Smooth linear statistics 

Theorem 3.1 (|34|) Let h be a C 2 -function on C with a compact support. 
Then 

Var Z r (h) = K + °^ \\Ah\\l, r -► oo , (3.2) 
where k is a positive numerical constant. The random variables 

r 2 

_ _ / 7,Jh\ - - 
y/k\\Ah\\ 2 

converge in distribution to the standard Gaussian law Af(0; 1) for r — > oo. 



r 2 f \ 
Z r (h) / h dm J 



Asymptotic formula (j3,2[) first appeared in Forrester and Honner jllj . 
It is worth mentioning that Theorem 13.11 persists for the other two CAZP 
models in the large intensity limit L — > oo |341 Part I]. 

It is instructive to compare (|3.2j) with the size of variations for a simple 
point process in the plane given by i.i.d. Gaussian perturbations of the 
lattice. Consider the point process 

S = {^(k + il) + rj k) i : (k, I) G I?} 

where r)k t i are independent standard complex Gaussian random variables. 
In this case, Var S r {h) ~ const ||V/i|||, for r — ► cxd. This is rather different 
from ()3.2|) . Asymptotic similarity to the flat CAZP Zf can be achieved by 
inventing special correlations between the perturbations % 4 In Section 0J 
we return to the idea of the fiat CAZP as a perturbed lattice. 
The proof of Theorem 13.11 starts with Green formula 



Z r (h)-— [ hdm = — [ \og\f* r \Ahdm, f*(z) 

IT J 2-7T J 



f(rz) 



y/Vaif(rz) 



The RHS is a non- linear functional on a Gaussian process /*. The rest is 
based on the method of moments a la Breuer and Major we expand 
the function Q ^ log \Q in Hermite polynomials in the space L^(e - '^ ) (the 
Wick expansion), and evaluate the moments of Z r {h) using the combinato- 
rial diagram technique. 

4 Lattice points are aggregated into clusters and each cluster scatters in a special 
(equiangular and equidistant) way |34l Part I, Introduction]. 
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3.2 Number of chaotic analytic zero points 



Let C C be a bounded domain with a piecewise smooth boundary. We 
are interested in the asymptotic behaviour of the random variable n(rQ) = 
Z r (tfi) for r — * oo. Forrester and Honner argued that the technique 
developed by Martin and Yalgin f22, for studying the Gibbs states of infinite 
systems of charged particles applied to the fiat CAZP gives 

Var n(rf2) = r ■ (r + o(l)) Length(<9r2) , r —* oo , 

r is a positive numerical constant. This is consistent with the idea that the 
variation of the number of points in ril should behave like the number of 
points in the 'strip' of constant size around the boundary d(rfl). 

Presumably, the method of Martin and Yalgin also yields that the ran- 
dom variables 

n(rO) — iT~ 1 r 2 m(Q) 
-y/r • r Length(<9r2) 

converge in distribution to 7^(0; 1) for r — ► oo. 

It would be interesting to find a counterpart of the law of the iterated 
logarithm; i.e. to find a function (j)(r) such that with probability one 

\n(r) — r 2 \ 
lim sup _ — = 1 . 

r-*oc V r 9y r ) 

Here n(r) = n ({\z\ < r}). 



3.3 The 'hole probability' and large deviations 

The next theorem proves an estimate conjectured by Yuval Peres: 
Theorem 3.3 ([34j) For r > 1, 

e" cir4 < P (n(r) = 0) < e~ c ^ , (3.4) 

where n(r) = n({\z\ < r}), and c\ and C2 are positive numerical constants. 

It would be interesting to check whether there exists the limit 

log" P (n(r) = 0) 
km 3 , 

r— »oo r 

and (if it does) to compute its value. 

The lower bound in (|3.4|) is obtained by an explicit construction. The 
upper bound follows from 
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Theorem 3.5 (|34|) For any 5 G (0, ?] and r > 1, 



n(r) 



1 



The proof of Theorem 13,51 uses tools from the entire functions theory. 
First, we show that with very high probability logmax r i} |/| is close to r 2 /2. 
Then, estimating log |/| from below, we show that with very high probability 
the average 

log|/(re^)|^ 



■2jt 



is also close to r 2 /2. From this, using Jensen's formula, we deduce Theo- 
rem I23J 

Theorems 13.31 and 13.51 are consistent with the results known for a one 
component Coulomb system of charged particles of one sign embedded into 
a uniform background of the opposite sign, Jancovici, Lebowitz and Magni- 
ficat [18;. It would be good to understand the asymptotic behaviour of the 
random variable r~ a [n(r) — r 2 ] for r — > oo and a > 1/2. At present, we 
understand the extreme case a = ^, and (not completely) the case a > 2. 
A plausible guess (motivated by JH]) is 



lim 



log logP (|n(r) 



> r a ) 



log r 



(2a -I, |<a<l; 
3a -2, 1 < a < 2; 
2a, a>2. 



In the first case, the normalized charge \n(r) — ^-| grows slower than the 
perimeter of the disc, in the second case it grows faster than the perimeter 
but slower than the area, in the last case (so called 'overcrowding') it grows 
faster than the area. According to the philosophy of ^H]) this should lead 
to a change of the asymptotic regime at a = 1 and a = 2. The technique 
we developed for the proof of Theorem 13.31 helps to analyze the case a > 2. 
The other two cases seem to require a different technique. 



4 Flat chaotic analytic zero points as a perturbed 
lattice 

How evenly do the flat CAZP spread over the plane? In the euclidean 
case we have a very natural system of points spread evenly throughout the 
plane, namely, the lattice points. Instead of comparing the random counting 
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measure n/ with its average -to, we consider the flat CAZP as a perturbed 
lattice {->fTi{k + il) + ^ j: k, I G Z}, for some (dependent) complex-valued 
random variables £k,i- We may hope for fast decay of the tail probabilities 
P(|£fc/| > A) f° r large A, uniformly in (k,l) £ Z 2 . The uniformity becomes 
trivial if the distribution of (£k,l) is invariant under lattice shifts. We treat 
the random variables ; as measurable functions on the space n = c z2 of 
two-dimensional arrays £: Z 2 — > C. 

Theorem 4.1 (|34|) There exists a probability measure ¥ on il, such that 

(i) P is invariant under the shifts of 7?; 

(ii) the random set S = {y/n(k + il) +£k,l : (k,l) E Z 2 } is distributed like 
the flat CAZP Z; 

(Hi) E ^e 6 ^ ' ' 2 ^ < cxd for some e > 0. 

This result gives no information about correlation between ^ Proba- 
bly, they can be chosen to be nearly independent on large distances. 

The proof of Theorem 14.11 does not use the Gaussian nature of the flat 
CAZP but only the uniform boundedness of the exponential moment of the 
'random potential' u(z) = log|/(z)| — \\z\ 2 . The main ingredients of the 
proof are the M. Hall's 'marriage lemma' (needed to match the flat CAZP 
with the lattice \/7rZ 2 ), and a potential theory lemma which can be useful 
in other discrepancy problems: 

Lemma 4.2 Let u be a bounded delta- subharmonic function on C (i.e., a 
difference of two subharmonic functions), and let Au = /i — m in the dis- 
tributional sense, ^ is a non-negative measure. Then for any bounded Borel 
set E C C 

H(E) < m(E +t ) and m(E) < fi(E +t ) , 

1/2 

where E+t = {z: dist(z, E) < t} is a t-vicinity of E, t = const ||u||oo • 

The boundedness of u is too strong for applications. It can be easily 
weakened by convolving u with a smooth convolutor supported by an ap- 
propriate disc. 

The two ingredients described above alone help to prove only a local 
result in the spirit of Theorem 14.11 Globalization is still a problem: after 
smoothing, the random potential u(z) = log|/(z)| — ^\z\ 2 is almost surely 
unbounded. Rare fluctuations appear somewhere on the infinite plane C, 
though probably far from the origin. To achieve some locality, we introduce 
a special adaptive metric p on C using Hdrmander's construction (141 Sec- 
tion 1.4]. This metric p is small where the potential is large. Then we use a 
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version of Lemma 14.21 making use of p-neighbourhoods instead of euclidean 
ones. 

It is worth mentioning that in Ajtai, Komlos and Tusnady studied 
high probability matchings of a system of iV 2 independent random points 
A = {Ai, ... X N 2} uniformly distributed in the square [0, N} 2 C R 2 with the 
grid {u>i, ...u N i} = [0, N) n Z 2 . They considered the average transportation 



the minimum is taken over all permutations tt of {1,2, ... N 2 }. Then with 
high probability 



For related results see Leighton and Shor ^3], and Talagrand [SSJ. In our 
global set-up we deal with infinite measures in the plane. Then, according 
to the lower bound in (|4.3|) . for any matching the average transportation 
distance tends to infinity in the N — * oo limit. This leaves no hope for the 
finite average distance matching between the Poissonian point process in R 2 
and a lattice, even without a quantitative estimate (iii). 

It would be interesting to find a hyperbolic counterpart of Theorem 14. II 



The hyperbolic CAZP, in contrast to the other models, depends on two 
intrinsic parameters: the mean density of zeroes per unit hyperbolic area, 
and the size of the sets on which we count the number of points. For in- 
stance, nL{D r ) is a number of the hyperbolic CAZP with intensity L in 
the hyperbolic disc of radius r. This leads to different asymptotic regimes, 
and makes the hyperbolic model richer than the other two. It is also worth 
mentioning that, for different values of the intensity L, the Gaussian ana- 
lytic function (|2.3j) exhibits completely different patterns in the asymptotic 
behaviour when z approaches the boundary of the unit disc Kahane |16l 
Chapter 13]. 

An interesting observation by Diaconis and Evans |S], and Peres and 
Virag [22] says that the real part of the hyperbolic random function (j2.3|) 
up to a constant term is a Poisson integral of a Gaussian random noise on 
the unit circle. In the case L = 1 this is classical white noise |29| p. 11], 
the case L = 2 corresponds to the Gaussian process on the Dirichlet space 




Ki<N 2 



const y/log N < T(A) < Const y^og N . 



(4.3) 



5 Hyperbolic CAZP 



flJ /2 (T) [HI Example 5.6]. 
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5.1 An exactly solvable model 

Here, we discuss a recent finding of Peres and Virag which pertains to 
the case L = 1. Recall that the /c-point correlation function p(zi, ... Zk) of a 
random point process is 

/ s r Pe(zi, ...Z k ) 

p{zi, ... z k ) = hm — — -sTjr — , 

where p e (zi, —Zk) is the probability that each disc {\z — zA < e}, 1 < 
j < k, contains at least one point of the process. For the random zero 
processes of a Gaussian analytic function, the limit on the RHS always 
exists. An equivalent definition says that p(z%, is a density of the 

/c-point correlation measure (|1.4|) with respect to the Lebesgue measure. 
The correlation measure also contains singular terms supported by the large 
diagonal; these terms are expressed via j-point correlation functions with 
j < k. Thus, the random zero process can be described by its correlation 
functions. 

Using Hannay's formulas Peres and Virag proved 



Theorem 5.1 (Peres- Virag |29j) The correlation function of the hyper- 
bolic CAZP with intensity L = 1 is 



p(z 1 , ... z k ) = tt n det 



1 



(1 - ZiZjf 



l<i,j<k 



This remarkable identity makes the hyperbolic CAZP with L = 1 an 
'exactly solvable model' among all CAZP. 5 In particular, it yields amazingly 
simple explicit expressions for the distribution of the number of zeroes n(p) 
in the disc {\z\ < p} and for the asymptotics of the 'hole probability': 

Corollary 5.2 (|29j) Let Z be the hyperbolic CAZP with intensity L = 1. 
Then 

(i) n(p) has the same distribution as YlJLi^j where {Xj} is a sequence of 
independent Bernoulli random variables with f(Xj = 1) = p 2 i ; 

(ii) for p — > 1 

P(n(p) =0) = exp I n 



(Hi) the ratio 



1-p 

n(p) — ~En(p) 
^Var n(p) 



(5.3) 



3 Peres and Virag observe that this is the only determinantal process among CAZP. 
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(with En(p) = j^-y and Varn(p) = -^n) converges in distribution to the 
standard Gaussian law J\f(0; 1) for p — > 1. 

In this case, Varn(/o) has the same order of magnitude as En(p) whilst 
in the flat case the variance grows only as a square root of the mean. This 
naturally reflects the difference between the hyperbolic and euclidean ge- 
ometries: in hyperbolic geometry the perimeter grows like the area, and 
much more random zeroes are located near the boundary circumference. 

We are not aware of counterparts of (ii) and (iii) for the hyperbolic CAZP 
with L^l. 

Loose ends 

In this lecture we've only touched the 'ground level' of the theory. There 
are plenty of interesting and deep developments. Among them are 

• scaling limits of zeroes of random polynomials (Ibragimov and Zeitouni 
|17j and Shiffman and Zelditch |31| ) and of random holomorphic sec- 
tions of high powers of Hermitian line bundles (Bleher, Shiffman and 
Zelditch @j); 

• solutions of random systems of algebraic equations, including sparse 
systems (see Edelman and Kostlan 0, Malajovich and Rojas |21j . 
Shiffman and Zelditch [32] and references therein); 

• distribution of real zeroes of random real polynomials (Maslova |2.3| . 
Dembo, Poonen, Shao and Zeitouni [5], Bleher and Di [3], Aldous and 
Fyodorov 0); 

• links with the zero distribution of chaotic eigenfunctions (Nonnen- 
macher and Voros Hannay |13| ) and with the distribution of 
eigenvalues of random matrices with independent complex Gaussian 
entries (Forrester and Honner , Diaconis and Evans [Hj , Dennis and 
Hannay and references therein), 

and each of them deserves a special lecture (cf. [SHI)- But these are different 
stories to be told by other people. 
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